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Abstract 



(^i Light- matter systems allow to realize a strongly correlated phase where photons are present. In these systems strong 
O^l correlations are achieved by optical nonlinearities, which appear due to the coupling of photons to atomic-like structures. 
$-h This leads to intriguing effects, such as the quantum phase transition from the Mott to the superfluid phase. Here, we 
address the two-dimensional Jaynes-Cummings lattice model. We evaluate the boundary of the quantum phase transition 
and study polaritonic properties. In order to be able to characterize polaritons, we investigate the spectral properties 
of both photons as well as two-level excitations. Based on this information we introduce polariton quasiparticles as 
appropriate wavevector, band index, and filling dependent superpositions of photons and two-level excitations. Finally, 
we analyze the contributions of the individual constituents to the polariton quasiparticles. 
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1. Introduction 

In the last few years proposals for new experimen- 
tal realizations of strongly-correlated many body systems 
emerged. Among them are ultracold gases of atoms 
trapped in optical lattices [H,[3;ll] and light- matter systems 
0, la S @l • The latter consist of light modes which are 
confined in coupled cavity arrays. Due to the finite overlap 
of their quantum mechanical wave functions, photons are 
able to tunnel between adjacent cavities. Free photons are 
non-interacting, however, an effective repulsive interaction 
can be achieved by coupling photons to atoms or atomic- 
like structures leading to interesting physical phenomena, 
which are subject to strong correlations. Different theoret- 
ical models have been proposed to describe these effects. In 
one scheme the atomic like structures are modelled by two 
level-systems, which leads to an interaction of the Jaynes- 
Cummings (JC) type 0,0. Another scheme is based on 
Jour-level systems [4[ and takes advantage of electromag- 
netically induced transparency Here, we focus on the 
first mentioned type of interaction. The physics of a single 
cavity at site i is described by the JC Hamiltonian 0] 



H( c = uj c a] a 



e at a', 



9(a>i °1 



a] aj 



), (1) 



where u> c is the resonance frequency of the cavity, i. e., the 
frequency of the confined photons, e is the energy spac- 
ing of the two-level system, and g is the atom-field cou- 
pling constant. The coupling between the atom and the 
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Figure 1: Schematic illustration of the JCL model in two dimensions. 
Red wavy arrows indicate photons, which are confined in optical 
cavities and interact with two-level systems (blue bubbles). 



photons is achieved by dipole interactions. The operator 
a\ (a i ) creates (annihilates) a photon at cavity i and at 
(a^ ) is the raising (lowering) operator of the two- level sys- 
tem. The JC Hamiltonian conserves the particle number 
hi = a] a- + at aj , which is a result of the rotating wave 



approximation 
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The full system of coupled cavities is 
modelled by the Jaynes-Cummings lattice (JCL) Hamilto- 
nian 



H 



JCL 



a] a„ 



(2) 



where the first term allows photons to tunnel between 
neighboring cavities i and j. The restriction to nearest 
neighbors is indicated by the angle brackets (• • • ) around 
the summation indices. The hopping strength t of the pho- 
tons is given by the overlap integral of the photonic wave 
functions. The second term in Eq. ([2]) describes the JC 
physics of each cavity and the last term with the chemical 
potential \i controls the total particle number N p = hi, 
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which is a conserved quantity of the JCL Hamiltonian. An 
illustration of the two-dimensional JCL model is shown in 
Fig-d Throughout this paper we use the dipole coupling g 
as unit of energy. Under these considerations, the physics 
of the JCL model depends on three distinct parameters: 
the hopping strength t, the detuning A = u c — e and the 
modified chemical potential (i — ui c . 

The fundamental excitations present in the JCL sys- 
tem are termed polaritons. Polaritons are superpositions 
of both photons as well as two-level excitations. At zero 
hopping strength t the JCL Hamiltonian reduces to the 
JC model with energies shifted by — ym% and can thus be 



solved analytically [3, [12J, From this solution it fol 



lows that the energy which is necessary to add two polari- 
tons to the cavity is larger than twice the energy which is 
necessary to add one polariton to the cavity, leading to a 
repulsive interaction between the particles. For small, non- 
zero hopping strength, the JCL system is in the so-called 
Mott phase, provided the polariton density is integer. In 
the case of light-matter systems the Mott phase can be re- 
garded as a state of frozen light. For t > two energies are 
competing: The kinetic energy which is gained by the hop- 
ping process of the photons and the repulsive interaction 
between the particles. For t larger than a critical hopping 
strength t c , it is energetically favorable for the particles 
to delocalize on the whole lattice and to Bose condense 
in a state of zero momentum. In this parameter regime 
the JCL system is in the superfluid phase. The critical 
hopping strength t c defines the boundary of the quantum 
phase transition from the Mott to the superfluid phase. 

Lately, there has been a great deal of research interest in 
understanding light-matter systems and in particular the 
JCL model. Most of the work has been devoted to study 
its ground state properties. The quantum phase transition 
from the Mott to the superfluid phase has been investi- 
gated by means of density matrix renormalization group 
13 , [3] j the variational cluster approach [TBI [l6| , quan- 
tum Monte Carlo (QMC) 17j, and analytically by means 
of strong coupling perturbation theory [l8|,[l|| . Results are 
also available on mean field level 0, [2(J ; some signatures of 
the quantum phase transition have also been determined 
from exact diagonalization of small systems consisting of 
a few cavities [g, [2l|, [22|. The spectral prop erties of the 
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JCL model have been evaluated in Rcfs. 

In the present paper, we study in detail the polaritonic 
properties of the JCL model. This work extends our recent 
results published in Ref. [l6| to two-dimensional systems. 
To be able to characterize polaritons, we first investigate 
the spectral properties of both particle species present in 
the JCL model, namely the photons as well as the two-level 
excitations. Based on this information we are able to intro- 
duce polariton quasiparticles as appropriate, wave vector, 
filling and band index dependent linear combinations of 
photons and two- level excitations. Furthermore we ana- 
lyze the contributions of the individual constituents to the 
polariton quasiparticles. For completeness we present the 
boundary of the quantum phase transition from the Mott 



to the superfluid phase as well and compare our results to 
QMC results obtained in Ref. [l7| . 

The remainder of this paper is organized as follows: in 
Sec. [5]we present the numerical method we are using. Sec- 
tion [3] contains the results for the quantum phase tran- 
sition and Sec. 2] provides a detailed investigation of the 
spectral properties of both particle species. Polaritonic 
properties are analyzed in Sec. [5] Finally, we conclude 
and summarize our findings in Sec. [6] 



2. Variational cluster approach 

We employ the variational cluster approach 0, [25| 
(VCA) as a numerical tool to study the quantum phase 
transition, spectral properties and polariton quasiparticles 
of the two-dimensional JCL model. In particular, VCA 
provides the single-particle Green's function G(k, ui) of 
the physical system H JCL and is based on the self-energy 
functional approach [26, 27 1 and the cluster perturbation 
theory 0, [29[ . VCA has been previously applied to light- 
matter systems in Refs. [Til, [l6|,|3(|. The main idea of VCA 



is that the self-energy Z of the physical system H JC is 
approximated by the one of a so-called reference system 
H' , which shares its interaction part with H JCL and is 
exactly solvable. Typically, H' is a cluster decomposi- 
tion of H JCL . The "optimal" reference system is deter- 
mined by requiring that an appro pria te functional fl[E] of 
the self-energy is stationary (see [26j for details). SI [E] is 
the grand potential of the physical system at the station- 
ary point. Additionally, at the stationary point of f2[E], 
Dyson's equation is recovered and thus the Green's func- 
tion of the physical system can be extracted there [26[ . To 
evaluate f2[E] numerically, the self-energy is parametrized 
by the single-particle parameters x of the reference system 
H' . Due to this parametrization the functional f2[E] be- 
comes a function f2(x), and is restricted to a smaller sub- 
space of self-energies. Physical quantities are evaluated 
at the stationary point of O(x) with respect to the single- 
particle parameters x. The accuracy of the results depends 
on the cluster size of the reference system, as correlations 
arc taken into account exactly on the cluster level. Thus 
convergence on physical quantities is achieved by increas- 
ing the cluster size of the reference system. 

The reference system W is solved by exact diagonaliza- 
tion, which we carry out by means of the band Lanczos 
method [3l|, |32[. To evaluate O(x) and the Green's func- 
tions G(k, oj) we apply the bosonic Q-matrix formalism 
[331 ]. In our previous work on the one-dimensional JCL 
model we adapted the VCA procedure such that it pro- 
vides the single-particle Green's function of both particle 
species, which involved as a subtlety the mapping of the 
two-level systems onto hard-core bosons. To summarize, 
we are able to extract the Green's function of photons 
G ph (\i, to) and the Green's function of two- level excitations 
G ex {k, to). From the Green's function G x (k, oj), where x 
stands either for photons (ph) or for two-level excitations 
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Figure 2: Phase boundary of the two-dimensional JCL model for 
zero detuning A = delimiting the first and the second Mott lobe 
from the superfluid phase. The phase boundary is evaluated for 
reference systems of size L. The inset shows the geometry of the 8 
site cluster. The marks refer to parameters where spectral functions 
arc evaluated. 



(ex), we obtain the single-particle spectral function 



A x (k, w) = ImG^k, w) 

7T 



(3) 



and the density of states 

N*(u)= /"^(k,a;)dk=i^^(k 1 a;). (4) 



3. Quantum phase transition 

The two-dimensional JCL model exhibits a quantum 
phase transition from the Mott phase, where polaritons are 
localized within the cavities, to the superfluid phase, where 
polaritons delocalize on the whole lattice 0, |T^, 17, ITc| . 
We evaluate the phase boundary for zero detuning A = 
by means of VCA from the minimal excitation energies of 
the system, see Fig.[2J Zero detuning means that the 
resonance frequency of the cavities is equal to the energy 
spacing of the two-level systems. We use the resonance 
frequency uj c and the energy spacing e of the two-level 
systems as variational parameters, i.e., x = {cj c , e}, see 
Ref. [l6| for details concerning the choice of the varia- 
tional parameters. Results converge quickly with increas- 
ing cluster size L. We evaluated the phase boundary for 
clusters consisting ofL = 2xl, L = 2x2 and L = 8 
sites, where the latter are arranged according to the inset 
of Fig. [21 In the case of L = 8 site clusters, the numerical 
determination of the stationary point of f2(x) turns out to 
be highly non-trivial. The reason for this behavior is sub- 
ject to further research, which is out of the scope of this 
work. The width of the Mott lobes decreases with increas- 
ing particle density, which results from the fact that the 
repulsive interaction itself depends on the particle density. 
The critical hopping strength £™, which specifies the tip 
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Figure 3: Photon spectral function A ph (k., lj), first row [(a), (e)], 
and density of states N ph (ui), second row [(b), ({)]. Two-level exci- 
tation spectral function A ex (k., ui), third row [(c), (g)], and density 
of states N ex (ui), fourth row [(d), (h)]. The spectral functions are 
evaluated for the parameters (a)-(d) t = 0.015, fi — u) c = —0.75, 
A = 0, marked with x and with the Roman numeral I in Fig.[2] and 
(e)-(h) t = 0.05, fi — ui c = —0.81, A = 0, marked with x and with 
the Roman numeral II in Fig. [2] Both parameter sets belong to the 
first Mott lobe. 
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of the nth Mott lobe, is t c w 0.055 for the first lobe and 
t\ « 0.0047 for the second lobe. The VCA result for t\ is in 
good agreement with the QMC data from Ref. [ItJ , where 
the authors obtained t\ ps 0.052 for the critical hopping 
strength of the first Mott lobe. 



4. Spectral properties 

Spectral functions A x (k, lo) and densities of states 
N x (uj) evaluated for both particle species are shown in 
Fig. [3] The spectral properties have been evaluated for 
the identical variational parameter set x = {uj c , e} as in 
the case of the quantum phase transition, for reference 
systems of size L = 2 x 2 and an artificial broadening 
0+ = 0.03. 

Spectral functions evaluated in the first Mott lobe con- 
sist of three bands. This can be understood from the zero 
hopping limit, where it is sufficient to consider a single 
cavity. The Hamiltonian then reduces to a block diagonal 
form, see for instance Refs. 0, |Tl| for a full solution of the 
single cavity problem. Each block corresponds to a spe- 
cific polariton number, i.e., the states contributing to one 
block are \n, \) and \n — 1, t)j where the first quantum 
number denotes the number of photons and the second 
one states whether the two-level system is excited (t) or 
in its ground state {\). Diagonalizing the block yields the 
dressed states \n, ±), where \n, — } is the ground state of 
the block and \n, +} its excited state. From these con- 
siderations it can be deduced that the single-particle ex- 
citations present in the spectral functions result from ex- 
citations from the ground state |n, — } to \n + 1, — ) or to 
|n + 1, +). This yields two distinct particle bands, which 
we denote as lu~ and uj^ , respectively. The same holds for 
the hole excitations, which corresponds to excitations from 
\n, — ) to | n — 1, — ) or \n — 1, +). We denote these bands 
as w^" and ut . However, for spectral properties evaluated 
in the first Mott lobe only one hole band exists, as the state 
|0, +) does not exist. This follows from the fact that in 
the zero hopping limit the block of zero particles consists 
only of one state |0, l). 

The location of the bands is identical for both the pho- 
tons and the two-level excitations, as they are connected 
by the dipole coupling. However, the weights differ signif- 
icantly. In particular, the particle band uj~ is much more 
pronounced in the photon spectral properties than in the 
two-level excitation spectral properties, whereas the oppo- 
site is true for the band which is almost not visible in 
A ph (k, ui) but carries significant weight in A ex (k, oj). The 
intensity of the hole band uj~^ is similar in both cases. 
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5. Polaritonic properties 

In this section, we investigate the polaritonic properties 
of the two-dimensional JCL model. To this end we follow 
Ref. fl6| , and introduce polariton quasiparticle and quasi- 
hole creation operators as appropriate superpositions of 
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Figure 4: Photon amplitudes /3 and two-level excitation amplitudes 
7 of the polariton quasiparticle creation operators (a), (c) p[_ . corre- 
sponding to the band uip and (b), (d) k corresponding to the band 

uip. In addition, the overlap A is shown in all panels (dashed line). 
The left column, (a)-(b), corresponds to the parameters t = 0.015, 
fi — oj c = —0.75, A = [see mark I in Fig. [2] and the right column, 
(c)-(d), to t = 0.05, fj,-LU c = -0.81, A = [see mark II in Fig.[2]. 



photons and two-level excitations 

pi,k = ^(k)4 + 7 p Q (k)a+, 
<k = ^(k)a k + 7 ?(k)<7k . 
The polariton quasiparticle states are defined as 



y / (^ob Q ,kPi jk IV , o) 

(ip Q I h aM hl k |-0o) 



and 



where polariton quasiparticle/quasihole operators are ap- 
plied on the iVp-particle ground state |^o}. In the follow- 
ing arguments we neglect for simplicity the band index a 
and wave vector k. The goal of this formulation is, to de- 
termine the coefficients of the linear combination j3 and 7 
such that the polariton quasiparticle states \ip p /h) describe 
best the true iV p ±l particle system. Therefore the overlap 
of \i> p /h) with the true (N p ± l)-particle eigenstates has to 
be maximized with respect to the coefficients 



max 

/3,7 



(f V±1 |^) 



This condition yields a generalized eigenvalue problem [If 
of the form 

A z = A S z , 
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where A is a 2 x 2 matrix containing the spectral weight of 
photons and two level excitations contributing to the spe- 
cific wave vector and band index, S is the overlap matrix 
and z determines the coefficients /3 and 7 up to a constant, 
which is specified by the conservation of the total spectral 
weight. The eigenvalue A determines the quality of the po- 
lariton picture as it corresponds to the normalized overlap 
specified above. Thus A is bound by [0, 1]. The polariton 
amplitudes j3 and 7 of the linear combination are shown 
in Fig. 0] for the particle bands uip. For the band uj~ the 
polaritons can be almost described by photonic excitations 
as P dominates over 7. In addition the polariton picture 
is well fullfilled as A « 1. For the band w+ the situation is 
reversed as 7 dominates over /?. In this case the polariton 
picture is solely modestly fullfilled as A w 0.3. Interest- 
ingly, for uj~ the amplitudes /3 and 7 are of opposite sign, 
whereas for cj+ they are of the same sign. This behav- 
ior is already observed in the limit of a single JC cavity 
and is a result of the special structure of the eigenvectors 
161 ] . However, as compared to the single-cavity limit, in 
the coupled-cavity system the modulus of the amplitudes 
is significantly altered. Finally, in both cases u>+ and uj~ , 
the polariton amplitudes /3 and 7 depend slightly on the 
wave vector k, which is also a pure lattice effect. 



6. Conclusions 

In the present paper, we investigated the polaritonic 
properties of the two-dimensional Jaynes-Cummings lat- 
tice model by means of the variational cluster approach. 
First, we evaluated the phase boundary delimiting the 
Mott from the supcrfluid phase and compared the criti- 
cal hopping strength which determines the position of the 
Mott lobe tip with quantum Monte Carlo calculations. 
We found good agreement between the two approaches. 
In addition, we evaluated spectral functions and densities 
of states for both particle species present in the Jaynes- 
Cummings lattice model, namely photons as well as two- 
level excitations. Spectral functions generally consist of 
four bands, two particle bands and two hole bands. An 
exception are spectral functions evaluated for parameters 
of the first Mott lobe, which contain two particle bands 
but only one hole band. Furthermore, we investigated in 
detail the difference in the spectral weight of photons and 
two-level excitations. Based on the spectral properties of 
both particle species, we were able to introduce polariton 
quasiparticles as wave vector, band index and filling depen- 
dent superpositions of photons and two-level excitations. 
The amplitudes of the linear combinations determine the 
character of the polaritonic quasiparticles. Moreover, the 
amplitudes depend on the wave- vector, which is a pure 
lattice effect. 
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